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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


UNIFORMLY  HIGH-ORDER  ACCURATE  NONOSCILLATORY  SCHEMES.  I* 

AMI  HARTENt  and  STANLEY  OSHERt 

Abstract.  We  begin  the  construction  and  the  analysis  of  nonoscillatory  shock  capturing  methods  for 
the  approximation  of  hyperbolic  conservation  laws.  These  schemes  share  many  desirable  properties  with 
total  variation  diminishing  schemes,  but  TVD  schemes  have  at  most  first-order  accuracy,  in  the  sense  of 
truncation  error,  at  extrema  of  the  solution,  (n  this  paper  we  construct  a  uniformly  second-order  approxima¬ 
tion,  which  is  nonoscillatory  in  the  sense  t'.iat  the  number  of  extrema  of  the  discrete  solution  is  not  increasing 
in  time.  This  is  achieved  via  a  nonoscillatory  piecewise-linear  reconstruction  of  the  solution  from  its  cell 
averages,  time  evolution  through  an  approximate  solution  of  the  resulting  initial  value  problem  and  an 
average  of  this  approximate  solution  over  each  cell. 


THE  DISCRETE  ONE-SIDED  LIPSCHITZ  CONDITION  FOR 
CONVEX  SCALAR  CONSERVATION  LAWS* 

YANN  BRENIERt  ano  STANLEY  OSHERt 

Abstract.  Physical  solutions  to  convex  scalar  conservation  laws  satisfy  a  one-sided  Upschiu  condition 
(OSLC)  that  enforces  both  the  entropy  condition  and  their  total  variation  boundedness.  Consistency  with 
this  condition  is  therefore  desirable  for  a  numerical  scheme  and  was  proved  for  both  the  Godunov  and  the 
Lax-Fricdrichs  scheme — also,  in  a  weakened  version,  for  the  Roe  scheme,  all  of  them  being  only  first  order 
accurate.  A  new,  fully  second  order  scheme  is  introduced  here,  which  is  consistent  with  the  OSLC.  The 
modified  equation  is  considered  and  shows  interesting  features.  Another  second  order  scheme  is  then 
considered  and  numerical  results  are  discussed. 


Some  New  Homogenization  Results  for  Discrete 
Boltzmann  Equations  with  Oscillatory  Data 

Thomas  Y.  Hou 

Department  of  Mathematics 
University  of  California  at  Los  Angeles 
Los  Angeles,  Ca  90024 


Abstract.  The  behaviors  of  multi-dimensional  discrete  Boltzmann  systems  with  highly 
oscillatory  data  are  studied.  Homogenized  equations  for  the  mean  solutions  have  been 
obtained.  Uniform  convergence  of  the  solutions  of  the  oscillatory  equations  to  the  solutions  of 
the  homogenized  equations  has  been  established.  Moreover,  we  find  that  the  weak  limits  of 
the  oscillatory  solutions  for  the  Broadwell  type  of  models  are  not  continuous  functions  of  the 
discrete  velocities.  Generalization  of  the  above  results  to  problems  with  multiple  scale  initial 
data  is  established. 


This  is  a  portion  of  the  author’s  Ph.D.  thesis.  Research  is  supported  by  ARO  Grant  No.  DAAG 
29-85-K-0190 
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Particle  Method  Approximation  of  Oscillatory 
Solutions  to  Hyperbolic  Differential  Equations 

Bjorn  Engquist  and  Thomas  Y.  Hou 


Department  of  Mathematics 
University  of  California,  Los  Angeles 
Los  Angeles,  CA  90024 


ABSTRACT 

Particle  methods  approximating  hyperbolic  partial  differential  equa¬ 
tions  with  oscillatory  solutions  are  studied.  Convergence  is  analyzed 
for  approximations  for  which  the  continuous  solution  is  not  well 
resolved  on  the  computational  grid.  Highly  oscillatory  solution  to 
the  Broadwell  and  variable  coefficient  Carleman  model  is  con¬ 
sidered.  Numerical  approximation  of  more  general  systems  are  dis¬ 
cussed.  Numerical  examples  are  given. 
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SOME  RESULTS  ON  UNIFORMLY  HIGH-ORDER  ACCURATE  ESSENTIALLY 
NONOSCILLATORY  SCHEMES 

Ami  HARTEN  * 

Department  of  Mathematics,  University  of  California,  Los  Angeles,  CA  90024,  U.S.A.,  and. 
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Department  of  Mathematics,  University  of  California,  Los  Angeles,  CA  90024,  U  S  A. 

Sukumar  R.  CHAKRA VARTHY 

Rockwell  Science  Center,  Thousand  Oaks,  CA,  US. A. 


We  continue  the  construction  and  the  analysis  of  essentially  nonoscillatory  shock  capturing  methods  for  the 
approximation  of  hyperbolic  conservation  laws.  These  schemes  share  many  desirable  properties  with  total 
variation  diminishing  schemes,  but  TVD  schemes  have  at  most  first-order  accuracy  in  the  sense  of  truncation 
error,  at  extrema  of  the  solution.  In  this  paper  we  construct  an  hierarchy  of  uniformly  high-order  accurate 
approximations  of  any  desired  order  of  accuracy  which  are  tailored  to  be  essentially  nonoscillatory.  This  means 
that,  for  piecewise  smooth  solutions,  the  variation  of  the  numerical  approximation  is  bounded  by  that  of  the  true 
solution  up  to  Q(h*-'),  for  0  <  R  and  h  sufficiently  small.  The  design  involves  an  essentially  nonoscillatory 
piecewise  polynomial  reconstruction  of  the  solution  from  its  cell  averages,  time  evolution  through  an  approxi¬ 
mate  solution  of  the  resulting  initial  value  problem,  and  averaging  of  this  approximate  solution  over  each  cell. 
To  solve  this  reconstruction  problem  we  use  a  new  interpolation  technique  that  when  applied  to  piecewise 
smooth  data  gives  high-order  accuracy  whenever  the  function  is  smooth  but  avoids  a  Gibbs  phenomenon  at 
discontinuities. 


1.  Introduction 

In  this  paper  we  consider  numerical  approximations  to  weak  solutions  of  the  hyperbolic  initial 
value  problem  (I VP) 

u,  +  /(“)*  =  0  •  +  a(u)ux,  (1.1a) 

u(x,  0)  =  u0(x).  (1.1b) 

Here  u  and  /  are  m  vectors,  and  a(u)  =  3//3«  is  the  Jacobian  matrix,  which  is  assumed  to 
have  only  real  eigenvalues  and  a  complete  set  of  linearly  independent  eigenvectors. 

The  initial  data  u0(x)  are  assumed  to  be  piecewise-smooth  functions  that  are  either  periodic 
or  of  compact  support. 

Let  v”  =  vh{Xj,  tn),  Xj  =  jh,  in  —  nr,  denote  a  numerical  approximation  in  conservation  form. 


*  Research  supported  by  NSF  Grant  No.  DMS85-03294,  ARO  Grant  No.  DAAG29-85-K.-0190,  NASA  Consortuim 
Agreement  No.  NCA2-IR390-403,  and  NASA  Langley  Grant  No.  NAGI-270. 
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Convergence  of  Vortex  Methods  for  Highly 
Oscillatory  Vorticity  Distribution 

Thomas  Y.  Hou 

Department  of  Mathematics 
University  of  California  at  Los  Angeles 
Los  Angeles,  Ca  90024 


Abstract.  Vortex  methods  for  two  dimensional  in  viscid,  incompressible  flow  with 
highly  oscillatory  vorticity  distribution  are  studied.  A  homogenization  result  for  Euler  equa¬ 
tions  is  obtained  up  to  a  finite  time.  Convergence  of  vortex  methods  in  the  case  when  continu¬ 
ous  vorticity  is  not  well  resolved  on  the  computational  grid  is  analylized.  Numerical  results 
are  given.  Comparisons  are  made  with  the  corresponding  finite  difference  approximation. 


This  is  a  portion  of  the  author’s  Ph.D.  thesis.  Research  is  supported  by  ARO  Grant  No.  DAAG 
29-85-K-0190 
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An  hierarchy  of  uniformly  high  order  accurate  essentially  non-oscillatory  shock, 
capturing  algorithms  was  developed.  Some  theory  and  numerical  experimentation  was 
done. 

A  correction  to  the  unsteady  full  potential  equation  for  flows  with  strong 
shocks  was  obtained.  This  modification  inputs  the  correct  entropy  jumps  at  shocks. 
Numerical  experiments  on  airfoils  were  sucessfully  performed. 

A  new  family  of  paraxial  wave  approximations  was  derived  and  was  applied  to 
computational  problems  in  seismology,  underwater  acoustics  and  artificial  boundaries. 
Theoretical  and  experimental  results  were  obtained.  The  family  also  included 
variants  of  parabolic  approximations  of  scalar  wave  equations. 

A  method  for  the  computation  of  highly  oscillatory  solutions  to  hyperbolic 
equations  was  obtained.  A  convergence  concept  which  makes  analysis  possible  in 
the  practical  situation  in  which  not  all  frequencies  are  well  resolved  is  developed. 
Convergence  of  an  average  approximation  is  established  for  a  general  class  of 
methods.  Applications  to  particle  methods  were  also  obtained. 
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Abstract 

An  entropy  correction  method  for  the 
unsteady  full  potential  equation  Is  presented. 
The  unsteady  potential  equation  Is  modified  to 
account  for  entropy  jumps  across  shock  waves. 
The  conservative  form  of  the  modified  equation 
is  solved  In  generalized  coordinates  using  an 
Implicit,  approximate  factorization  method.  A 
flux-biasing  differencing  method,  which 
generates  the  proper  amounts  of  artificial 
viscosity  In  supersonic  regions,  is  used  to 
discretize  the  flow  equations  in  space. 
Comparisons  between  the  present  method  and 
solutions  of  the  Euler  equations  and  between  the 
present  method  and  experimental  data  are 
presented.  .  The  comparisons  show  that  the 
present  method  more  accurately  models  solutions 
of  the  Euler  equations  and  experiment  than  does 
the  isentroplc  potential  formulation. 

Nomenclature 

speed  of  sound 
metrics  of  coordinate 
transformation.  Equation  (7) 
airfoil  chord 
time  step  In  computational 
space 

indices  of  grid  points 
Identity  matrix 
Jacobian  of  coordinate 
transformation 
reduced  frequency  based  on 
semi  chord 
Mach  number 
total  speed 

distance  from  point  of  rotation 
to  a  point  in  the  flow  field 
gas  constant 
entropy 
time 

contravarlarrt  velocities, 
Equation  (6} 
cartesian  coordinates 
angle  of  attack 

p2"y/J 

ratio  of  specific  heats 
Jump  In  potential  across  the 
wake 
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difference  operator 
computational  coordinate 
di recti ons 
polar  angle 
airfoil  pitch  angle 
density 

p/J 

biased  density 

computational  time 

velocity  potential 

average  of  quantity  across  the 

wake 


isentroplc 
normal  to  shock 
steady  airfoil  location 
pitch  axis  location 
poi nts  on  ai rf oi I 
minimum  airfoil  pitch  angle 
maximum  airfoil  pitch  angle 
free  stream  conditions 


lower  side  of  wake 
iteration  number 
upper  side  of  wake 
sonic  conditions 


Introducti on 

Modern  aircraft  typically  operate  at  high 
speeds  where  aeroelastlc  Instabilities  are  more 
likely  'to  occur.  To  successfully  predict  and 
analyze  such  aeroelastlc  phenomena,  the  aircraft 
designer  needs  methods  that  accurately  predict 
the  aerodynairKc  loads— steady  ado  unsteady— that 
the  vehicle  experiences.  Since  many  critical 
aeroelastlc  phenotnqna  occur  ay  transonic  speeds, 
methods  based  on  NJjnear  Aerodynamic  theory 
cannot  accurately  predict/  many  aeroelastlc 
responses.  Thus,  it  Necessary  to  use  an 
aerodynamic  method  that  /,  can  predict  time- 
accurate  solutions  of  nonlili^ar  flows  and  that 
can  accurately  model  shock  'tjaves  and  thei r 
unsteady  motions. 

/ 

t 

When  shock  waves  appear  in  transonic  flow 
fields,  aerodynamic  loads  predicted  using 
potential  flow  theory  can  be  grossly  inaccurate 
and  even  Aultlvalued.  Multiple  solutions  of  the 
potential  equation  were  first  observed  in  two 
dimensions  by  Steinhoff  and  Jameson. i  Salas 
and  Gumbert?  showed  that  the  phenomenon  is  not 


On  Numerical  Particle  Methods 


Yizhao  Hou 


Introduction 

Particle  Methods  have  received  a  great  deal  of  attention  in  the  past  few 
years.  In  Fluid  Dynamics,  Vortex  Methods  have  been  successfully  used  to  simu¬ 
late  High  Reynold  number  fluid  flows.  They  were  first  introduced  by  Rosenhead 
[33]  and  were  subsequently  developed  by  Chorin  ,  Leonard  among  other  contri¬ 
butors  (see  Rcf[25]  for  a  detailed  bibliography  ).  Applications  of  Vortex  Methods 
include  turbulent  combustion([l5],[34]),  aerodynamic  calculations([26],[9]),  flows 
through  heart  valves  [25]  and  flows  of  variable  density  [1]. 

In  Physics,  on  the  other  hand.  Particle  Methods  have  been  used  for  numeri¬ 
cal  solutions  of  kinetic  equations  such  as  Boltzmann  equations,  Vlasov  equations. 
They  are  currently  used  for  studying  a  number  of  physical  phenomena,  for 
instance  the  lasermatter  interaction  in  inertial  confinement  fusion  problems.  See 
the  book  of  Hockney  and  Eastwood  [18]  for  more  information  along  this  line. 

It  has  been  observed  in  practice  that  Particle  Methods  are  often  used  for  the 
problems  in  which  the  number  of  particles  is  too  small  for  the  traditional  asymp¬ 
totic  error  estimate  to  be  applicable,  and  yet  very  useful  results  are  still  possible 
(  see,  eg.  [30]  )  .  The  previous  mathematical  analysis  for  Particle  Methods  (  see, 
eg.  [16],  [4],  [5]  and  [32]  )  have  not  been  able  to  give  a  satisfactory  explanation 
to  this  feature  of  Particle  Methods  . 
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ABSTRACT 


Different  variants  of  parabolic  approximations  of  scalar  uave  equations 
are  derived  and  their  properties  are  analysed.  These  equations  are  of  general 
form  which  includes  the  ones  used  in  seismology,  underwater  acoustics  and 
other  applications.  A  new  version  of  the  parabolic  approximation  is  derived 
for  heterogeneous  media.  It  has  optimal  properties  with  respect  to  wave  reflec¬ 
tion  at  material  interfaces.  The  amplitudes  of  the  reflected  and  transmitted 
waves  depend  continuously  on  the  interface.  Existence,  uniqueness  and  energy 
estimates  are  proved. 


ABSTRACT 


A  new  family  of  paraxial  wave  equation  approximations  is  derived. 

These  approximations  are  of  higher  order  accuracy  than  the  parabolic  approxi¬ 
mation  and  they  can  be  applied  to  the  same  computational  problems  e.g  in 

seismology,  underwater  acoustics,  and  as  artificial  boundary  counditions.  The 

(€<3 

equations  are  written  as  systems  which  simplif/  computations.  The  support  and 
singular  support  are  studied,  energy  estimates  are  given  which  prove  the  well- 
posedness.  The  reflection  and  transmission  are  showed  to  be  continuously  depen 
ding  on  material  interfaces  in  heterogeneous  media. 


COMPUTATION  OF  OSCILLATORY  SOLUTIONS 
TO  HYPERBOLIC  DIFFERENTIAL  EQUATIONS 


BJORN  ENGQUIST1 


Abstract.  Difference  approximations  of  hyperbolic  partial  differential  equations 
are  studied.  The  error  is  analyzed  in  approximations  for  which  the  continuous 
solution  is  not  well  resolved  on  the  computational  grid.  A  convergence  concept  is 
introduced  with  the  purpose  of  making  analysis  possible  and  describing  the  practi¬ 
cal  situation  in  which  not  all  frequences  in  the  solution  are  well  resolved.  Con¬ 
vergence  of  an  averaged  approximation  for  linear  problems  is  established  for  a 
general  class  of  methods.  For  special  methods  and  a  weakly  nonlinear  problem  it 
is  showed  by  numerical  experiments  and  in  one  case  by  analysis  that  the  discrete 
approximation  converges  to  the  solution  of  the  homogenized  equation. 
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Efficient  Implementation  of  Essentially 
Non-Oscillatory  Shock  Capturing  Schemes 


Chi-Wang  Shu 
and 

Stanley  Osher 


In  the  computation  of  discontinuous  solutions  of  hyperbolic  conservation 
laws,  TVD  (total-variation-diminishing),  TVB  (total-variation-bounded)  and  the 
recently  developed  ENO  (essentially  non-oscillatory)  schemes  have  proven  to  be 
very  useful.  In  this  paper  two  improvements  are  discussed:  a  simple  TVD  Runge- 
Kutta  type  time  discretization,  and  an  ENO  construction  procedure  based  on 
fluxes  rather  than  on  cell  averages.  These  improvements  simplify  the  schemes 
considerably  —  especially  for  multi-dimensional  problems  or  problems  with 
forcing  terms.  Preliminary  numerical  results  are  also  given. 


Fronts  Propagating  with  Curvature  Dependent  Speed: 

Algorithms  based  on  Hamilton-Jacobi  Formulations 

Stanley  Osher 
and 

James  A.  Sethian 

We  devise  new  numerical  algorithms,  called  PSC  algorithms,  for  following 
fronts  propagating  with  curvature-dependent  speed.  The  speed  may  be  an 
arbitrary  function  of  curvature,  and  the  front  can  also  be  passively  advected 
by  an  underlying  flow.  These  algorithms  approximate  the  equations  of  motion, 
which  resemble  Hamilton-Jacobi  equations  with  parabolic  right-hand-sides,  by 
using  techniques  from  the  hyperbolic  conservation  laws.  Non-oscillatory 
schemes  of  various  orders  of  accuracy  are  used  to  solve  the  equations,  providing 
methods  that  accurately  capture  the  formation  of  sharp  gradients  and  cusps  in 
the  moving  fronts.  The  algorithms  handle  topological  merging  and  breaking 
naturally,  work  in  any  number  of  space  dimensions,  and  do  not  require  that  the 
moving  surface  be  written  as  a  function.  The  methods  can  be  also  used  for 
more  general  Hamilton-Jacobi- type  problems.  We  demonstrate  our  algorithms  by 
computing  the  solution  to  a  variety  of  surface  motion  problems. 
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Abstract 

In  the  computation  of  discontinuous  solutions  of  hyperbolic  conservation  laws, 
TVD  (total-variation-diminishing),  TVB  (total-variation-bounded)  and  the  recently 
developed  ENO  (essentially  non-oscillatory)  schemes  have  proven  to  be  very  useful. 
In  this  paper  two  improvements  are  discussed:  a  simple  TVD  Runge-Kutta  type 
time  discretization,  and  an  ENO  construction  procedure  based  on  duxes  rather 
than  on  cell  averages.  These  improvements  simplify  the  schemes  considerably  - 
especially  for  multi-dimensional  problems  or  problems  with  forcing  terms.  Prelimi¬ 
nary  numerical  results  are  also  given.4A_ 
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We  devise  new  numerical  algorithms,  called  PSC  algorithms,  for  following  fronts  propagating  with 
curvature-dependent  speed.  The  speed  may  be  an  arbitrary  function  of  curvature,  and  the  front 
can  also  be  passively  advected  by  an  underlying  flow.  These  algorithms  approximate  the  equations 
of  motion,  which  resemble  Hamilton- Jacobi  equations  with  parabolic  right-hand-sides,  by  using 
techniques  from  the  hyperbolic  conservation  laws.  Non-oscillatory  schemes  of  various  orders  of 
accuracy  are  used  to  solve  the  equations,  providing  methods  that  accurately  capture  the  forma¬ 
tion  of  sharp  gradients  and  crisps  in  the  moving  fronts.  The  algorithms  handle  topological  merging 
and  breaking  naturally,  work  in  any  number  of  space  dimensions,  and  do  not  require  that  the 
moving  surface  be  written  as  a  function.  The  methods  can  be  also  used  for  more  general 
Hamilton- Jacobi-type  problems.  We  demonstrate  our  algorithms  by  computing  the  solution  to  a 
variety  of  surface  motion  problems. 
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Uniformly  High  Order  Accurate  Essentially 
Non-oscillatory  Schemes,  III 

Ami  Harten 
Bjorn  Engquist 
Stanley  Osher 
Sukumar  R.  Chakravarthy 


We  continue  the  construction  and  the  analysis  of  essentially  non-oscillatory 
shock  capturing  methods  for  the  approximation  of  hyperbolic  conservation  laws. 

We  present  an  hierarchy  of  uniformly  high-order  accurate  schemes  which  generalizes 
Godunov's  scheme  and  its  second-order  accurate  MUSCL  extension  to  an  arbitrary 
order  of  accuracy.  The  design  involves  an  essentially  non-oscillatory  piecewise 
polynomial  reconstruction  of  the  solution  from  its  cell  averages,  time 
evolution  through  an  approximate  solution  of  the  resulting  initial  value  problem, 
and  averaging  of  this  approximate  solution  over  each  cell.  The  reconstruction 
algorithm  is  derived  from  a  new  interpolation  technique  that,  when  applied  to 
piecewise  smooth  data,  gives  high-order  accuracy  whenever  the  function  is  smooth 
but  avoids  a  Gibbs  phenomenon  at  discontinuities.  Unlike  standard  finite 
difference  methods  this  procedure  uses  an  adaptive  stencil  of  grid  points  and, 
consequently,  the  resulting  schemes  are  highly  non-linear. 
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We  continue  the  construction  and  the  analysis  of  essentially  non-oscillatory  shock  capturing 
methods  for  the  approximation  of  hyperbolic  conservation  laws.  W'e  present  an  hierarchy  of 
uniformly  high-order  accurate  schemes  which  generalizes  Godunov  s  scheme  and  its  second- 
order  accurate  MUSCL  extension  to  an  arbitrary  order  of  accuracy  The  design  involves  an 
essentially  non-oscillatory  piecewise  polynomial  reconstruction  of  the  solution  from  its  cell 
averages,  time  evolution  through  an  approximate  solution  of  the  resulting  initial  value 
problem,  and  averaging  of  this  approximate  solution  over  each  cell.  The  reconstruction 
algorithm  is  derived  from  a  new  interpolation  technique  that,  when  applied  to  piecewise 
smooth  data,  gives  high-order  accuracy  whenever  the  function  is  smooth  but  avoids  a  Gibbs 
phenomenon  at  discontinuities.  Unlike  standard  finite  difference  methods  this  procedure  uses 
an  adaptive  stencil  bf  grid  points  and.  consequently,  the  resulting  schemes  are  highly  non¬ 
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